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Abstract 

To describe structural peculiarities in inhomogeneous media caused by the 
tendency to the close packing of atoms a formalism based on the using of the 
Riemann geometry methods (which were successfully applied lately to the 
description of structures of quasicrystals and glasses) is developed. Basing on 
this formalism we find in particular the criterion of stability of precipitates 
of the Frank-Kasper phases in metallic systems. The nature of the "rhenium 
effect" in W-Re alloys is discussed. 

Recently the interest grows to the significance of crystallogeometrical factors, first of all 
connected with the close packing of atoms, for the structure of condensed matter [lj . There is 
a whole class of metallic systems which structure is apparently determined by the tendency 
to the close packing of ionic spheres with different radii including quasicrystals, metallic 
glasses 0, small metallic clusters || and Frank-Kasper (FK) phases ||. It is worthwhile to 
stress that the considerations of the closest local packing may contradict with the existence 
of long-range crystal order. It distinguishes essentially the real three-dimensional case from 
the two-dimensional one where (for the case of spheres of equal radii) the triangle lattice 
provides the closest packing both for the whole space and its every part i.e. both globally and 
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locally [|J. Classical problem about the distribution of spheres in three-dimensional space 
to obtain the closest packing has not yet solved rigorously (it is the part of the 18th Gilbert 
problem). However, it is known that spheres may be packed locally with higher density 
than in fee, hep and other close packed lattices [Q. Therefore the situation of geometric 
frustrations arises when the optimum structure from the point of view of local surrounding 
cannot be optimum globally. The concept of frustrations appeared to be useful in particular 
for the consideration of the structure of disordered systems 

The cause of the existence of the geometric frustrations is the fact that the Euclidean 
space cannot be occupied completely by regular tetrahedra (the regular tetrahedron is the 
structure unit providing the closest packing of four atom groups). An elegant technique has 
been proposed in Ref.6 to construct the structures of FK phases when the latter are obtained 
from the closest packing of tetrahedrons in Riemannian space by the introducing of the net 
of structure disclinations (SD) decurved the space and then the filling of the Euclidean space 
by the structure units arisen. Here we will use a similar approach to analyze the problem 
of formation of heterogeneous state (HS) in metallic alloys. Such state is the subject of 
great interest of researchers in the field of material science lately. From the experimental 
point of view it is distinguished by the peculiar diffuse X-ray scattering corresponding in the 
real space to the presence of relatively small clusters containing only hundreds or even tens 
of atoms JF|. The examples are the so called athermic cj-phase in some Ti- and Zr-based 
alloys || and in Cr^-jAla, 0, a-phase in Fei- x Cr x |H|, precipitates of W 3 Re phase with 



A15 structure in Wi-^Re^ |[L1| . Recent investigations of Ti^a-Fe^ alloys by Moessbauer 



effect [12| have shown the presence of geometric frustrations i.e. the difference of the type 
of short-range order in local surrounding of Fe nuclei from the long-range order in a crystal 
as a whole. The hypothesis has been proposed in Ref. 13 about the existence of small 
icosahedral clusters in this alloys. 

Generally speaking, the structure of crystal phases is determined by the whole number of 
factors. Among them, apart from the considerations of close packing, the spatial orientation 
of covalent bonds, peculiarities of electronic structure near the Fermi surface etc are the 



most frequently discussed. Therefore it is naturally to expect that the nature of HS may be 
different for different systems. In particular Krivoglaz [[!]] stressed the role of peculiar features 
of the shape of the Fermi surface for the formation of HS in Ti- and Zr-based alloys. This 
mechanism is however hardly to be universal since in the whole number of systems with HS 
the electron mean free path is small and therefore the alloy smearing of the Fermi surface 
has to be essential. 

The alternative mechanism is proposed in the present work basing on the considerations 
of the closest packing of atoms for small atomic groups and geometric frustrations connected 
with them which may lead to the formation of HS with characteristic scale of order of few 
interatomic distances. The model proposed is rather rough and pretends to describe only 
some main features of the phenomenon. Nevertheless it stress, to our opinion, a peculiar 
role of the factor which is important for many real systems but was not taken into account 
in previous works. 

Consider the precipitates of A15 phase in bcc host (e.g. of W 3 Re phase in the solid 
solution Wi_ x Re x [11]) as an example of heterogeneous state of the type under consideration 
because A15 structure is the simplest and most well-known representative of FK phases [|TJ. 
It contains eight atoms per cell, two of them having icosahedral coordination and six of 
them being surrounded by polyhedra with 14 vertices || . The latter may be obtained from 
the regular icosahedron by the introducing of the edge disclination with the axis passing 
its center and the Frank vector -2tt/5. Using such a procedure it was shown in Ref.6 that 
the net of SD for A15 structure consists from three mutually orthogonal sets of equidistant 
parallel disclinations with the period being equal to the lattice constant a. Other FK phases 
differs from A15 structure only by the geometry of SD net. 

To solve the question about the stability of HS one needs to generalize the approach . 
For this aim we use the basic relation between the curvature tensor Rijki and the disclination 
density tensor B\j pA 

Rijkl — ~ e ijp e klq^q , (1) 
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where e^ p is the unit antisymmetric tensor. In light of Eq. (1) the space curvature providing 
the ideal tetrahedral packing in the approach is created by the introducing of the partial 
edge disclinations with Frank vector Q = 7°20' which compensate the deficiency of dihedral 
angles at the packing of five tetrahedra around the common edge ||16|| . To construct FK 
phases it is necessary to "decurve" the space by the introducing of SD with the average 
density 9^ D = —9[f where the line means the average over the volume containing large 
enough number of cells. Thus, in contrast with the approach of Ref.6 where SD were 
introduced in curved space we introduce two sets of disclinations in Euclidean space with 
the mean total density being equal to zero. At the same time local variation of the curvative 
of the lattice (i.e. the deviation of local disclination density from zero) is possible and 
connected with the elastic distortions of the bonds. The corresponding stresses are similar 
to that which have discussed in Ref.5 for glasses. 

Apart from the energy of elastic distortions the electronic ( "chemical" ) contribution to 
the energy describing the tendency to the closest packing of atoms also exist. If in some place 
the density of SD differs from the mean density in the FK phase and therefore according to 
(1) the lattice appears to be locally curved the local values of the packing density will be 
higher or lower than in the FK phase. 

In the framework of the approach considered the energy of the FK phase counted from 
the energy of the "host" phase with 6 = may be written as 

E = J dr/[%(r)]+^, (2) 

where / is the density of "chemical" energy depending on the degree of the atomic packing 
which is connected with the disclination density tensor 9ij, E e i is the elastic energy. As usual 
]I7J the latter may be represented in linear elasticity theory in the following form 

Eei = \j J drdr'rf»{v)H vpKT {v - r')r) KT {r'), (3) 

where H vpKT is the Green tensor for internal stresses, r\ vp is the incompatibility tensor de- 
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scribing the density of sources of internal stresses. According to Ki\ it is denned as 



V ^ 2 ppcpLvq,p @vp ^upq^pq,p ® pv) i (4) 

where a vp is the disclination density tensor. Separating the singular part of the Green tensor 
we may represent the energy E e i as the sum E el = Ef + E^ of the energy of distortions 
in the disclination cores E { ° ] and the energy of elastic deformations outside the cores E el . 
When dislocations are absent one has r\ vp = — l/2{6 up + pu ). For single disclination with 



the Frank vector (0,0,f2) the tensor r\ vp = —6 up has the only non zero component [14 



= -Q5(p), (5) 

where p = (x,y). In the case of the only set of edge disclinations considered here the tensor 
8ij according to (5) may be characterized by the only scalar parameter 9 = Tr8ij which 
equals to the mean value of the Frank vector in a given point. It is proportional to the 
scalar curvature Rijij (see (1)). For the given geometry the energy of disclination cores (per 
unit length) may be represented as 



E. 



(o) 

el „ / j2 „n2 



L 



e d 2 P 6 2 (p) 



where eo is the energy of the core of the disclination with unit Frank vector determining 
by the relation e S(p) = ^H 3333 (0). Note that for the elastic continuum e = and it is 



necessary to use the quasicontinuum model 17 to describe the core energy correctly. 



It is a common practice in the material science to treat inhomogeneous state mainly as 
the multiphase state with the coexistence of regions with different crystal structures each 
of them may be stable in principle in the whole space. In this case the inhomogeneity may 
be described in terms of elastic distortions due to the conjugation of crystal lattices on the 
boundary of the precipitate and the host. To calculate the distortions the Eshelby model [I5[ 



is usually used, their energy appearing to be proportional to the volume of the precipitate. 
It is the consequence of the fact that in the Eshelby model the deformation e\j connected 
with the phase transition is constant inside the precipitate. Then, according to Eshelby, 



the deformation e\j may be simulated by a system of disclination loops on the boundary 
of the precipitate. Since the energy of every loop is proportional to its radius R (equal to 
precipitate radius) and the number of loops is proportional to R 2 their total energy turns 
out to be proportional to R 3 . In our model the precipitate has the type of short-range order 
(e.g. icosahedral) which cannot correspond to long-range order for any bulk crystal phase. 
Therefore in contrast with the Eshelby model the deformation appears to be inhomogeneous 
not only outside the precipitate but also inside it. As the result it will be shown below 
that the elastic energy of the precipitate grows with the increase of R faster than R 3 and 
therefore the stabilization of HS becomes possible. 

The real form of the function f{6) in (2) is unknown. In accordance with our choice of the 
zero point for the energy we have for host phase (9 = 0) /(0) = and the minimum of f(9) 
lies at the value 9 = 9q corresponding to the closest tetrahedral packing in the Riemannian 
space [f|. We suppose for the simplicity 

f = a(9 - 29 )9, (0 <fl<# ), (6) 

where the parameter a > depends on the explicit form of the interatomic interactions. 
Generally speaking, the function / may contain also terms proportional to (V9) 2 , but in 
the framework of the variational approach used below (see (8)) they do not influence on the 
results. 

To demonstrate the possibility of the formation of heterogeneous state we use the direct 
variational approach and construct explicitly the distribution 9(r) ^ const leading to the 
lower value of the total energy than the homogeneous state. The problem may be solved in 
the simplest way for the cylindrical precipitate of the close packed phase. Let Oz be the axis 
and R the radius of the precipitate, the axis of all the disclinations are parallel to Oz and 
9 = 9(p) where p = \/x 2 + y 2 . 

Since it is obvious that the discontinuities on the boundary between the host and pre- 
cipitate are energetically unfavorable we restrict ourselves by the case of their coherent con- 
jugation when the discontinuities are absent. In this case the disclinations cannot abrupt at 
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the boundary and should form loops with external segments in the host. In contrast with 
the homogeneous A15 phase the "polarized" distribution of disclinations is typical for our 
case namely the uniform distribution of positive partial disclinations inside the precipitate 
and the "cloud" of negative disclinations outside it. From topological considerations (the 
vanishing of the total Frank vector) one has 

/ d 2 P 6(p) = (7) 

To evaluate the minimum of the energy (2) we use the trial function 



0(p) 



0i, p<R 

Sgggsr, R<p<R + A (8) 
0, p>R + A 



answering the requirement (7) automatically. The value A describes the thickness of the 
"cloud" and should be found from the solution of the variational problem. In the model 
under consideration the energy of elastic distortions concentrated in the disclination cores 
in the cloud would play the role of the surface energy. However, due to the condition (7) 
the total number of disclinations in the cloud is equal to that in the precipitate itself and is 
proportional to its volume (in 2D case ~ R 2 ). Therefore in contrast with the Eshelby model 
it is impossible to separate explicitly the bulk and surface parts from the precipitate energy. 
We use for H 3333 (p — p') = H(p — p') the known expression |l] 

2p p 2 



H(r) 



l-ln-P- 

Rr 



+ C (9) 



1 - uStt 

corresponding to the continual approximation in the elastically isotropic model where p is 
the shear modulus, v is the Poisson ratio, R c is the cutoff radius of order of the size of a 
crystal, C is the constant which does not contribute to the energy E^' for a given geometry 
(total Frank vector equals zero). To describe the energy of disclination cores E^ correctly 
one should use the quasicontinuum approximation JT7]] 
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where kd is the Debye wave vector, 7 is the parameter depending on the explicit form of 
phonon dispersion curves in the model. 

Then the elastic energy (3) per unit length of the precipitate does not depend on the 
cutoff radius R c and has the form 

L =el T^^ R > A) (10) 

where 

jtoA) = Ri{A + R)2 , 

192a 4 A 2 (A + 2R) 

A (5 A 3 + 20 A 2 R + 26 A R 2 + 12 R 3 ) 

A + Rv 
R ). 

The minimizing of the function $(R, A) with respect to A gives the dependence of the cloud 



-12 R 2 (A + R) 2 In (- 



thickness A = A(R). As a result the total energy may be represented in the following form 

f = irR 2 ap(6 l - 29 )9 1 + 6f^—^(R) (11) 
Lj 1 — v 

where ip(R) = ip(R, A(R)), 9 = o /p,p = 1 + 2e /a. The minimum of energy (11) corre- 
sponds to the disclination density 

9l = 1 + *P(R)/kR 2 ' (12) 
where the precipitate radius R is found from the equation 

Jl^=2(2*(fl)+«fl») (13) 

Here 

5f(l - u)np 

and the value Sf = a9 2 describes the chemical energy gain (per unit volume) at the transition 
to the ideal tetrahedral packing. 

As it follows from Eq. (13) the precipitate radius R depends on the parameter k char- 
acterizing the ratio of the chemical contribution to the elastic one. According to (14) its 
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value varies from 0.1 to 1 at Sf as 0.01 — 0.1 eV per atom. Direct calculations show that 
in all these limits of the variation of n the function if)(R) is almost constant at R < 2a and 
increases sharply at R > 3a. On that ground, the equilibrium value of R grows slowly with 
k increase (see Fig.l). At k < 0.1 the continual approach under consideration does not hold 
because R < a. In this range of values of k the disclination density 6 X determined by Eq. 
(12) is less than the minimal possible value corresponding one disclination with Frank vector 
7°20' per precipitate. 

The spatial distribution of the dilatation obtained by numerical integration is shown in 
Fig. 2. One may see that the internal part of the precipitate is expanded and the dilatation 
en vanishes sharply at the boundary. The mean value in the precipitate < En > increares 
approximately linear with the precipitate radius. This demonstrates the qualitative differ- 
ence between the model under consideration and the Eshelby model. Since the average 
deformation increases with R increasing the elastic energy will grow faster than R 2 in 2D 
case ( or than R 3 for spherical precipitate). 

Putting the equilibrium values of &i,R in Eq.(ll) one obtains the expressions for the 
energy change at the formation of the precipitate 

E = s f ,#h = - wL^m <0 (15) 

L J O #0 - 0i 1 - v 
Thus in the model under consideration the precipitates are stable and the heterogeneous 
state has the lower energy than the bcc lattice. Accordingly to the dependence 9 1 on R 
its value (12) , generally speaking, is not coincide with one of the values corresponding 
homogeneous FK phases. 

Thus we show that for certain ratio between "chemical" and elastic properties of the 
media (« > 0.1) the heterogeneous state may be energetically favorable owing to the crys- 
tallogeometrical factors (tendency to the closest local packing). It is worthwhile to stress 
again the most important features of this state in the framework of the model under con- 
sideration (i) sizes of the precipitates are of order of few lattice constants (ii) the structure 
of precipitates is distorted in comparison with the structures of FK phases in an infinite 
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crystal (iii) precipitates are surrounded by the cloud of the distorted host lattice with the 
thickness of order of the sizes of the precipitates. The existence of this cloud has an obvious 
sense: since it is impossible to fill the whole space by the regular tetrahedra an attempt to 
realize such filling locally in some part of the space results inevitably to the appearance of 
voids at the boundary. The introducing of the external segments of disclination loops lead 
to the elimination of these voids and coherent conjugation of the precipitate with the host 
surrounded. 

In the specific case when the host phase is bcc a sharp increase of the solubility of light 
interstitial impurities is one of the characteristic features of the mechanism considered of 
the formation of HS. It is connected with the presence of large tetrahedral voids in close 
packed structures and their absence in bcc lattice. This feature together with the change of 
the type of short-range order and the appearance of elastic distortions can help to identify 
this kind of HS experimantally. 

The results obtained here may explain some important features of the so called "rhenium 
effect" (the improvement of the ductility of W and Mo at the doping by Re). According to 
Ref. 11 the increase of the solubility of interstitial impurities connected with the appearance 
of the precipitates of the FK phases in particular W 3 Re is the key feature in this phenomenon. 
The results presented here allows us, on the one hand, to understand the causes of the 
appearance of the precipitates and, on the other hand, to point out the specific mechanism 
of their influence on the solubility of interstitial impurities. Indeed, considerable dilatation 
in the cloud surrounded the precipitate may lead to the essential energy gain at the transfer 
of the interstitial atom from the host to the shell. It may prevent the appearance of the 
carbide precipitates or impurity segregation at the intergrain boundaries. 

Note in the conclusion that the transition from the two-dimensional case considered 
here to the three-dimensional one is not trivial since to "dress" the spherical precipitate 
by the cloud of compensated defects we cannot restrict ourselves by the disclinations only 
and it is necessary to introduce the density of disclination loops as an additional variable. 
However the results presented here seem to be sufficient to demonstrate the instability of 
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the homogeneous state under some conditions. 

To resume it is worthwhile to note that the present work show the interconnection of some 
specific facts known in the material science with non-trivial properties of three-dimensional 
space and therefore demonstrate the general physical meaning of these facts. 

This work is supported by Russian Basic Research Foundation, grants 95-02-05656 and 
95-02-06426. We are grateful to A.M.Kosevich for helpful discussions. 
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FIGURES 

FIG. 1. The dependence of the equilibrium values of precipitate radius R, disclination density 
9\ (in units of 6q) and the value of ip in Eq.(ll) on the k parameter (14). 

FIG. 2. The distribution of the dilatation e^k in the percipitate (solid line) on the polar 
coordinate r (R=4, A = 1), and the dependence of average dilatation (efcfc) on the precipitate 
radius R (dashed line) in the units #i(l — 2i/)/(27r(l — v)). 
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